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We consider the entropy and decoherence in fermionic quantum systems. By making a Gaussian 
Ansatz for the density operator of a coUection of fermions we study statistical 2-point correlators and 
express the entropy of a system fermion in terms of these correlators. In a simple case when a set of 
A'^ thermalised environmental fermionic oscillators interacts bi-linearly with the system fermion we 
can study its time dependent entropy, which also represents a quantitative measure for decoherence 
and classicalization. We then consider a relativistic fermionic quantum field theory and take a mass 
mixing term as a simple model for the Yukawa interaction. It turns out that even in this Gaussian 
approximation, the fermionic system decoheres quite effectively, such that in a large coupling and 
high temperature regime the system field approaches the temperature of the environmental fields. 

PACS numbers: 



I. INTRODUCTION 

The density operator contains complete information about quantum statistical systems, and hence it can be used 
to study various properties of such systems, such as correlators, particle numbers, entropy and decoherence. However, 
the evolution of realistic physical systems is governed by interacting field theories, and only rarely the density operator 
is known beyond a perturbative approximation, which can be, nevertheless, very useful for weakly coupled regimes. 

Even the Gaussian part of the density operator contains important information about the entropy and decoherence of 
the system, which can be neatly encoded in the statistical two point function, as was firstly pointed out independently 
by two groups of authors [U This correlator approach to entropy, decoherence and classicalization has been 
extensively used in the context of weakly interacting bosonic systems [TUIO] . However very little is known about the 
entropy and decoherence in fermionic systems, and the corresponding literature is scarce fllfll3|. In this paper we 
present a first study of entropy and decoherence in relativistic fermionic field theories. For simplicity, we consider here 
only simple bilinear interactions, which are in field theory known as mass mixing. Since our Hamiltonian is quadratic in 
the fields, an initial Gaussian density operator will remain Gaussian as the system evolves, and a complete information 
about the density operator can be given in terms of equal-time 2-point correlators, which is the strategy we use in 
this work. For pedagogical reasons, we begin by considering coupled fermionic quantum oscillators, and only then 
move on to field theory. 

If the density operator of a system p is known, the (information) entropy S'vn can be calculated by the von Neumann 
formula, 

= -{Hp)) = -Tr[/51n(p)] . (1) 

Now, by making use of the Heisenberg evolution equation for the density operator, one can easily show that the 
von Neumann entropy is conserved for closed systems. In practice however no observer O will have access to a 
complete information of any nontrivial system S (with many interacting degrees of freedom), making the system 
open. Such systems will interact with an environment E which is, by definition, inaccessible to O. The loss of 
information associated with this inaccessibility generically leads to decoherence jT4Ul8] , a rather qualitative concept 
that describes how a system evolves into a state which most closely resembles a classical state. However, at the 
same time the (reduced) von Neumann entropy ([l]) of the system alone ^ is no longer conserved due to this loss of 
information. Entropy generation thus provides a quantitative measure of decoherence and classicalization. 
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^ I'he von Neumann entropy of the reduced system density operator is in literature also known as the entanglement entropy |19l 1201 . as 
the information about the entanglement between the system and environment is lost in the reduced density matrix. 
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This statement requires clarification, because decoherence is an observer-dependent concept, whereas entropy ([T]) 
can be defined without introducing an observer, but of course it does depend on the system - environment split. To 
make our work as general as possible, we kept the observer implicit throughout. Yet, a natural observer O is the 
one which entangles with those states that diagonalise the density operator. In this work we identify those states 
as the Fock states for which statistical (average) particle number is defined. Decoherence is now induced by tracing 
over the inaccessible environmental degrees of freedom, and is perfect from the viewpoint of the observer entangled 
with the Fock states that does not see the entanglement between system and environment. As these Fock states get 
more occupied, and the system's entropy increases, the system gets more classical. This observer plays a special role 
in the class of all observers sensitive to Gaussian properties of the density operator, and that is that this observer 
sees the most quantum properties of the system. All other observers, such as the momentum or position operator, 
will perceive the system as more classical. We can make this more concrete by considering the example of a highly 
squeezed pure state (that describes e.g. linearized cosmological perturbations at the end of inflation). That state will 
be perceived as classical by the observer that measures the position spread of the state (in the sense that Aa; will be 
much greater than in the pure vacuum state), but it will be perceived as a pure quantum state by the observer that 
entangles with the Fock states that diagonalize the system's density operator. 

But, what can be the inaccessible information that can be justifiably called an environment? The most common 
example is a thermal bath of particles interacting weakly with the system that is observed. In this case, the system- 
environment correlations ~ also known as entanglement - are not observable. These correlations can be in 2-point SE 
correlations (such as considered in this work) or in higher order n-point functions. Higher order correlators are always 
suppressed by some power of the coupling constant, and hence they are typically small in a weakly coupled regime. 
For example, in Refs. 019] the leading SE correlator that is neglected corresponds to a 3-point SE correlator. 

In this work we take the interaction to be of a Yukawa type, schematically 

-Cint ~yiji't(j)-ipj , (2) 

where (j) and denote a scalar and fermionic quantum field, respectively. The simplest approximation in which one 
can treat this interaction is to neglect quantum fluctuations of the scalar field, i.e. to replace the scalar field by its 
expectation value, 

Within this Gaussian approximation the Yukawa interaction reduces to a mass mixing term. 

Ant -'ipiMijipj , Mij = yij(j), (3) 

the Hamiltonian becomes quadratic in the fields, and the problem becomes exactly soluble. We shall use numerical 
techniques to obtain an exact solution to this simplified problem, and we shall express the Gaussian density operator 
in terms of equal-time 2-point correlators. The information we consider inaccessible to O is in the SE and EE 2-point 
correlators, and hence entropy gets generated. The so-called Gaussian von Neumann entropy for the system field, 
which is derived from the Gaussian density matrix alone, yields a good quantitative measure of decoherence for nearly 
Gaussian systems. However, for highly non-Gaussian systems, one has to modify the entropy definition to incorporate 
the relevant non-Gaussian features of the state j8j. In the simplified problem Q the Gaussian von Neumann entropy 
can be analytically calculated in terms of 2-point (statistical) correlators of the system degrees of freedom. This has 
also been done for various bosonic systems in Refs. [3 H]. In this work we set out to derive the Gaussian (von 
Neumann) entropy for fermionic systems in terms of correlators of the system. 

The Gaussian entropy is generally not conserved in the presence of interactions, which could be either environ- 
mental interactions, or self-interactions. Several case studies for bosonic systems [T] [21 [71 |H1 [lO] have indeed shown 
that the Gaussian entropy increases for interacting systems, thereby quantitatively describing decoherence and clas- 
sicalization. As far as we know, a quantitative description of decoherence for fermionic systems is still lacking. A 
better understanding of decoherence in fermionic systems can be applicable in many situations of physical interest. 
For instance, one species of fermions could mix with others through mass-type terms, such as quarks through the 
CKM matrix [21] [52] or neutrinos through the PMNS [231 IMj matrix. Other examples include Yukawa interactions, 
or condensed matter systems with interacting fermions. Here a framework is provided for calculating the growth 
of entropy for a system of fermions interacting via a fermionic mass matrix ([s]). As explained above, this model 
represents the simplest (Gaussian) approximation to the more realistic Yukawa interaction 

The outline of this work is as follows: in section [H] the simplest example of a one-dimensional fermionic harmonic 
oscillator is discussed. A general Ansatz is made for the density operator, after which the particle number and entropy 
are derived in terms of the statistical correlators of the system. In section [lI A| we make a connection with some existing 



literature by working with the density operator in the coherent state basis. Next, in section II B the simplest possible 
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interactions are added to the fermionic system: N environmental fermionic oscillators coupled bilinearly to the system 
oscillator. Though not completely realistic, this example provides an insight into how a loss of information leads to 
an increase in the entropy of a fermionic system, and some specific examples are shown. In section |III| we switch 
our attention to the more realistic fermionic quantum field theory. After discussing diagonalisation of the Dirac 
Hamiltonian in III A an Ansatz is made for the density operator in terms of mixing particle and antiparticle states 
in section [ill B| The Gaussian entropy is derived in terms of the statistical correlators. In section III C| thc fermionic 
entropy is generalised in the presence of N fermionic degrees of freedom. Finally, in section [HID the first realistic 



example of entropy generation in fermionic quantum field theory is discussed, which is the simple case of one fermionic 
species mixing with other species through mass terms. 



II. ENTROPY GENERATION IN FERMIONIC QUANTUM MECHANICS 



The most general Ansatz for the density operator of a free fermionic quantum mechanical system (fermionic harmonic 
oscillator) with the Lagrangian, 

Lo ^ ^\idt - Lo{t))^ , (4) 

can be written as ^, 

p{t) = exp(-ai/>^?/>) , (5) 

where a{t) is a (complex valued) function of time, and \/Z is the normalisation constant determined by the usual 
trace condition, 

Tr[p(t)] - 1 , (6) 

and "0 is the (Grassmannian) fermionic operator (here expressed in the Schrodinger picture) satisfying the usual 
canonical anticommutation relation, 

{^A^} = 1. (7) 

Now making use of the Grassmannian nature of the operators ip'^ and ip with (-0)^ = 0, (V"^)^ = 0, and of Q, we can 
expand ([5| as, 

p(<) = l(l+[c-'^-l]iv) , (8) 

where we introduced the fermionic number operator N = -ip^tp with TV" — N {n — 1,2, ..). The Hilbert space of this 
theory is two dimensional, and can be conveniently represented in terms of the Fock space basis vectors {|0), |1)}, 
defined by, 

N\n) ^ n\n) . (9) 
The trace of the density operator Q is easily evaluated in this basis, 

such that the general Gaussian fermionic density operator is properly normalised according to Q by Z = l + exp(— a). 
It is also convenient to express the density operator as 

m = - n{t)) + (Mt) - l)N , (11) 



^ The anticommutation relation ^ implies that the other possible Gaussian term exp{—b^ip^ ) can be expressed in terms of exp(—aip^ ip) 
plus an appropriate change in the normalisation constant, and hence does not constitute a new term. In presence of interactions, the 
lagrangian can be written as = Lq + ^inti where Ljnt = — J^V' — ip^j^ji- In this case the density operator can still be written as 
in where now ■4> and denote the suitably shifted fields, as shown in appendix B 
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where the average particle number n is defined as, 

{N) = Tr[p{t)N] 

The (von Neumann) entropy is then simply, 

-Tr[pHp)] 



nit). 



(12) 



S 



{n\ I [(1 - n) + {2n - 1)N] In [(1 - n) + (2n - 1)N] | \n) 



n=0,l 



This evaluates to 



S 



-(1 — n) ln(l ~ n) — nln(n) . 



(13) 



which is the standard expression for the entropy of n free (non-interacting) fermions, where h is the average number 
of fermions in the system defined in (12 1. For an analogous discussion of a bosonic oscillator we refer to appendix [A| 
Let us now make a connection with the familiar expressions for a thermal fermionic density matrix [25) . According 
to the Fermi-Dirac distribution, the average occupancy of a state with energy E is given by, 



1 



(14) 



where /3 = l/iksT) is the inverse temperature. Of course, 

npD = (N) = Tr[pth^] 



where pth denotes a thermal density operator. By comparing with the general expression for p (11) and with (12), it 
is now easily seen that the thermal density matrix is obtained upon identification, a — >■ /3E, such that, 



Pth 



,I3E 



1 - c 



fiE 



N 



or, equivalently, 



Pth = (1 - nth) + (2nth - . 



(15) 



(16) 



which is, as expected, of the same form as the general entropy (13). The thermal density operator (161 implies the 
following well known expression for the entropy of a thermal Fermi gas, 

•S'th = -(1 - "th) ln(l - nth) - nth In(nth) • (17) 



Let us now consider a bit more closely Eq. ( 13 ). In the spirit of the Schwinger-Keldysh out-of-equilibrium formalism, 
it is convenient to introduce the following 2-point functions, 

ts++{t;t')= {Timi^^m 
is+-{t-t') = -{^\t')m) 

iS'+{t-t')^ (^(t)^nO) 

is-{t-t')^ {fm)^\t')]), 



(18) 



where here ■0(0 ri^) &de the Heisenberg picture operators, and T and T denote time ordering and anti-time 
ordering operations, defined as. 



iS++{t; t') = e{t - t')iS-+{t; t') + e{t' - t)iS+-{t; t') 
iS'^it; t') = e{t - t')iS+-{t] t') + 9{t' - t)iS-+{t- 1') , 



(19) 



such that. 



iS 



++ 



■iS- 
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The retarded and advanced Green functions are then, 

iS' = iS++ - iS+- = -{iS— - tS-+) , iS" = tS+- 
Notice that tS'^ and iS^ can be also written as 

The statistical and causal (spectral) two point functions arc defined as 



eit'-mmJHt')})- 



(20) 



1 



1 



1 



(i,t'))-^([^(t)>nO]) 



(21) 



such that p^(t;t) = 1/2 {S''{t;t) = i). By making use of the identity, = {l/2)['tp\'ijj\ + (l/2){?/>1', -0}, one can 
obtain a simple relation between n{t) and F^{t] t): 



nit) = = ^-F4t;t) 



1- A^(t) 



where in the last step we defined 



A^{t) = 2F^{t-t) = 1 - 2n{t) = tanh(^) . 



(22) 



(23) 



Note that in Eq. (22 1 the average particle number is computed using time dependent operators ?/'(t), ■0^(i), although 
n in Eq. (12) was computed using Schrodinger picture operators a-nd the density operator ([5]). The time 

dependence of the operators can be absorbed under the trace into the density operator, which then takes the form 
([5]) for equal-time operators. 

Eq. ( 22 ) represents a relation between the invariant of the correlators and the Gaussian invariant of the density 
matrix, which are in this single fermion case simply F^(t;t) and a{t), respectively. In different systems with multiple 
correlators and a more complicated Gaussian density matrix such a relation can still be found. An example is the 
bosonic case, discussed in footnote [3] and in appendix [A} or the fermionic field theoretical case, discussed in section 

imi 



Note that n e [0, 1] and e [—1, 1], which can be appreciated from Eq. (23 1, making the interpretation of as 
the invariant phase space area of the state for fermions dubious. It is hence better to think about A^ as the Gaussian 
invariant of a fermionic state, while n = (1 — A^)/2 g [0, 1] is more like the phase space area. 

For thermal states, for which hth € [0,1/2], A^th acquires natural values, A^th G [0,1], and hence there is no 
problem. In fact, A^ becomes negative only when higher energy states are overpopulated, i.e. when n > 1/2. 
Relation (23 1 allows us to relate the fermionic entropy (13 1 to the Gaussian invariant A^(i), 



1 + A^, / 1 + A ^, \ 1-A 



In 



(24) 



which is to be compared with the analogous expression for bosons in Eq. (AlO) of appendix [A| 



A. Coherent states 



In order to make a connection to the existing literature |11H131E5] . here we rephrase our results in terms of fermionic 
coherent states 16*), defined by, 

m = 9\e) . (25) 



^ The definition | |23^ is the fermionic equivalent of the invariant (phase space area) A of a bosonic Gaussian state in Eq. | |A8[ | of 
appendix|A] where we present an analogous derivation of the entropy for bosons. The form of the Gaussian invariant for fermions 
is so simple because the fermionic density operator is diagonal in TV = ^/'t-i/i, implying that the fermionic density matrix is diagonal in 
the fermionic particle state basis. 
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When expressed in terms of Fock (particle number) states the coherent ket and bra states are given by, 

\e) = \o)-e\i), (0| = (o|-(i|^", (26) 

where 9 — 9* and Grassmann variables obey a Grassmann algebra, 9i6j = ~9j9i (recall that complex conjugation for 
Grassmann variables is reminiscent of a hermitian conjugation, {9i9j)* = 9*9* ~ —9*9*). By making use of the well 
known relations, 

^|0)=0, ^^|0) = |1), V'|1> = |0), V^^|i)=o, 



one sees that (26 1 is indeed an eigenstate of the operator ip with the eigenvalue 9. Note that the Fock space element 



|0) commutes with Grassmann variables, while |1) = "0^10) anticommutes, such that the coherent states \9) commute 



with Grassmann variables. The coherent state \9) is fixed uniquely by the requirement (25) up to a normalisation 



constant = 1 + b99, where 5 is a complex number. Our choice of normalisation corresponds to 

{9\9) = l + 99, (27) 
which is Grassmann valued. One may attempt to normalise to unity by choosing = 1 — 99/2. The problem with this 



is that then the operator ■0^ does not act on \9) in a desired manner. In fact, one can show that the normalisation (27) 
is uniquely fixed by the requirement, 

i'^ie)^-^^. (28) 



Indeed, when acts on 16*) defined in (26 1 one gets and when the derivative —d/d9 acts on the same state, one 
again gets A different normalisation would not give this result. Finally, as a final check of consistency we consider 
how the anticommutator acts on \9), 

^n\o) - ( - ^1 + ^^O) 19) ^ {^^^ - .^t) I,) . (iL, + e±) 19) = 19) , 

as it should be from ([t]). When projected on a coherent state basis, the elements of the density matrix (111 become 
of the form, 

p{9\9;t) = {9'\p{t)\9) = {l-n)+n9'9 = {l-n)exp(^-^-^9'9^ , (29) 

which is not diagonal. This is to be contrasted with a diagonal Ansatz used e.g. in Ref. flSl. 

Let us now consider properties of the coherent state basis in more detail. Taking a trace of the density operator in 
the coherent state representation yields,^ 



Tr [p] ^ J d-9 J d9exp{99){9\p\9) = J ^9 J d9{l + 99) [(1 - n) + n99] = 1 , (30) 
where in the last step we used the usual integration rules, J d9 — 0, J d99 — 1. The integration measure factor exp{99) 



in ( 30 ) is necessary to get the traces correctly. 
One can now use a decomposition of unity ^, 



le^ d9 d9exp{-99)\9){9\, (31) 



The expression for the trace in Eq. l |30| can be derived as follows: the trace of an operator O is in the Fock basis defined as 
TrP] = (OlOlO) + (llOll), where the Fock space elements can be expressed in terms of coherent states as 



|o> = J dee\e), \i) = - j de\e), (o| = j dee{e\, (i| = j 



^ Note that Tr[/e] = 2, as it should be, where Ig is given in Eq. |3l]l. 
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to recast p as, 

p= I d9'd9'e-'^''^'\e'){e'\p j dOdOe 

dO'de'c-""' I d0d9e-'^''p{e',9,t) 



d9'd9'e^'^''^'\9') [(1 -n)+ n9'9] j d9dee-'^\9\ 



(32) 

where p{6' , 9,t) are elements of the density operator in the coherent state representation (see Eq. (p9|), 

p{0',e;t) = \e') [{l-n)+ne'9] {9\ - \9')p{9\9;t){9\ , (33) 



with p{9', 9; t) = Z^^ ey.^{M9'9) given in Eq. p{9', 9; t) of Eq. ^ is obviously non-diagonal. However p of Eq. 

(34) 



(32) can be cast in a diagonal basis by inserting the Ansatz 

p{e,e';t)= f dCrfCe'"«+fV(c) , 



such that 



P= / dCrfCIC>nC)(C|. 



(35) 



This is the so-called Glauber P representation [57] for fermions. Inverting (34) the function P{C) is related to the 
density matrix in the diagonal elements of the coherent state basis as 



F(() = / d9d9e-'^^-^^ {9\p\i 



(36) 



The elements of the density matrix in the coherent state basis have been found in (29), and by integrating over 9,9 
in (36) one finds: 



P(C) = n+(l-n)CC. 



(37) 



It is possible to return to the Fock basis via the P representation (35) using Eqs. (26) and (37) and integrating over 



m = E - ") + (2" - = |0)(1 - n){0\ + |l)n(l| 



(38) 



n=0 



For one degree of freedom the (diagonal) Fock number basis is by far superior to the coherent state basis ( 35 ) for 
studying properties of the fermionic density operator, an important example being the von Neumann entropy defined 
in ([T]) and calculated in (13). The reason is that the Fock states are orthogonal, contrary to the coherent states which 
satisfy {C\9) = e^^ , see Eq. (27). Still, the von Neumann entropy can be derived from the density operator in the 
coherent state basis by using the replica trick, which is demonstrated in appen dix [P} There we also generalise to N 
fermionic degrees of freedom, a case which is discussed in more detail in section |III C[ 



B. Fermionic interactions in quantum mechanics 

Interactions can be included in the quantum mechanical fermionic theory ^ by introducing general current terms 
into the Lagrangian 

= Lo+ Lint ; ^int = - V'^ j'V' • (39) 

Formally the linear current terms can be absorbed into the free field theory ^ by shifting the fermionic fields, see 
appendix [B] For these free shifted fermionic fields the von Neumann entropy is conserved. 

In realistic situations it is very hard to have a complete information about the current operator jjp{t) however. 



which makes the diagonalisation procedure (B7) impracticable, or even impossible. Namely, in condensed matter 
systems, the coupling current is often given by a superposition of many (fermionic) degrees of freedom, whose precise 
time evolution is not known. In a quantum field theoretic setting one can have for example a Yukawa coupling term. 
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—y4>(x)ip(x)'4'{x), such that th e co upUng current corresponds to a composite operator, j^{x) = y4>{x)'ip{x), making 
the diagonahsation procedure (B7) very hard, if not impossible. For that reason we adopt here the point of view that 
no (useful) information is known about the evolution of the current . This loss of information leads to entropy 
generation, which is what we study next. 

The simplest nontrivial example is the quantum mechanical case when the current consists of N environmental 
oscillators in thermal equilibrium. In this case. 



N 



(40) 



where the ipq. represent the environmental fermionic oscillators. The form of the current (40) is motivated by mass 
mixing, which can be considered as an approximation to the Yukawa coupling, cf. Eqs. (|2H3| 



The system is 

represented by a single fermionic oscillator ^pJ: which is coupled bilinearly to the environmental oscillators through 
couplings Xi. The interaction between the environmental oscillators is assumed to be zero in our toy model. The loss 
of information in this case is that we cannot observe (correlations of) the environmental oscillators, nor its interaction 
with the system. The complete action of system, environment and interactions in our toy model is 



(41) 



with 



N 
i=l 

JV 



(42) 



Note that by the hermiticity of Lint, all A* = \i are real. The fermionic oscillators only depend on time, i.e. 'ipx = i'xit) 
and tpq. — ^g. (t) . The anticommutation relations satisfied by ipx and ipq. are 

{Vi,,(t),V'Ut)}-%, 



with all others being zero. Of our interest are the statistical correlators (21 1, which are for our model defined as 





-\< 


mt)M{t')]) 


F,^,,{t;t') 


-\'^ 


[^,,w>;(i')]) 




-\'^ 




^qiX (^1 ^ , 




[%At)Mt')])^ 



(43) 

where il)x{t) and tpqi{t) are here Heisenberg picture operators. Our goal is to calculate the entropy for the system. For 
a free fermionic theory the entropy is given by Eqs. (23) and (p4|). Without interactions, the Gaussian invariant A is 



constant and the entropy is conserved. If we switch on interactions the Gaussian invariant and the entropy change in 
time. The entropy of the system with the Gaussian Ansatz ([5| for p with time dependent a{t) is related to Axx as in 
Eq. E^: 



Sx{t) 



l + Axx{t), fl + Axx{t)\ 1-Axxit), (l-Axx{t) 



In 



In 



(44) 



with 



Axx{t) = 2Fxxit;t) = 1 - 2nxxit), 



(45) 
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and fixxit) the average particle number for the system fermions. The proper way to derive the entropy of the system is 
to trace over the environmental degrees of freedom in the density operator, and calculate the entropy from this reduced 
density operator. The corresponding reduced von Neumann entropy is the same as Eq. (44), i.e. S'^^(t) = Sx{t).^ 
Thus, in order to investigate the growth of entropy for the system fermionic oscillator ipx, we should find the statistical 
correlator Fxx(t,t) defined in Eq. (43 1. The equations of motion for the fermionic operators follow from the action 

(Ell), 

N 





- UJo)lpx{t) 


1=1 






= ^dx{t) 






N 
i=l 









(47) 

From these equations of motion we can derive coupled differential equations for the statistical equal-time correlators 

N 

Y^X, {Fq^x{t;t)-Fxq^it;t)) 



idtFxx{t\t) 
{idt - (wo - ^j))Fxqj {t; t) 

{idt - (ujj - ujo))Fqjx{t; t) 
iidt - {uj,~ujj))Fq^g^{t;t) 



N 



-^]Fxx(t; t) + ^ KFq.g^ (t; t) 



N 



\jFxx{t]t) -Y,^^Fq,<lMt) 



i=l 



XiFxq- {t; t) — XjFq.x{t; t) 



(48) 



These conditions can be solved with suitable initial conditions. We take the system fermionic oscillator to be initially 
in a state with average particle number zero. The environmental fermionic oscillators are assumed to be in a thermal 
state according to the Fermi-Dirac distribution with energy Ei = uJi, see Eq. (14). Thus 



Fxx{t0]tQ) — - 

Fq^q. (to; to) ^^ij^ tanh ^ 
Fxqi [to] to) — Fq.x{to; to) = . 



2 



(49) 



With these initial conditions Eqs. (48) can be solved. We first treat the simple case of two coupled oscillators, then 
the general case of N coupled oscillators. 



Two coupled fermionic oscillators 



For two coupled fermionic oscillators we consider the case = 1 in the action (41 ). Thus there is one environmental 
oscillator ijjq = 7/;^^ coupled to the system oscillator ipx through a coupling A = Ai. This simple example can be solved 



® The proof goes as follows. The reduced density matrix is defined as p'^'^'^ = TrE[p], where the subscript E denotes the environment, 
which in this example is the group of oscillators {V'<ji }■ The reduced von Neumann entropy is the usual 5^'f^(t) = -Tr[p^°'^ lnp'^<=<^]. Now 
most importantly, if we calculate correlators of the system, we have 

(1 - A,,)/2 = (i/;t^4 = Tr[p^t^,] = J2 K|KJ..(n,„|pVii^.;|n<,„>-|n9i>|n4 =Trs[/3'°'^Vit^^]. (46) 

Here lipx, ipt and p are taken to be in the Schrodinger picture, which ensures that ij)x'4'x does not evolve with time, whereas p controls the 
evolution. The trace is taken over a complete set of orthogonal, time independent Fock states. Eq. ( |46| shows that the correlators, and 
thus the invariant area | |45| l and entropy \AA\ for the system are the same whether you first trace over the environment in the density 
matrix or you consider the full density matrix. 
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Figure 1: System entropy as a function of ujot for N = 1 
environmental oscillator. The parameters are uji = 1.5a;o, 
A = O.Sojo and P = [luq)^^ . The entropy oscillates be- 
tween zero and the thermal entropy (17l (dashed line). 
The maximum entropy Smax ~ lii(2) is indicated with the 
dotted line. 



Figure 2: System entropy as a function of uiot for N = 1 
environmental oscillator. The parameters are lji — l.Sojo, 
A = O.SoJo and /3 = 0.1(ujo)~^ ■ In this case the en- 



tropy almost reaches the thermal entropy (17 1 (dashed 
line), which is in turn almost equal to the maximum value 
•Smaa; = ln(2) (dotted Imc) . 



analytically. The procedure is explained in appendixjC] As a final result we find an explicit expression for the Gaussian 
invariant of the system represented by the fermions ip^, 



^xx{t) = l-2nE 



where 



2A 



sm 



-{t-to) 



(50) 



He 



V (^0 



4A2 



(51) 



with the frequencies wq and ioi of the system and environment oscillators, respectively. The entropy of the system 
Sx is subsequently found using Eq. ( 44 1 . Figures [T] and [2] show the evolution of entropy for a system coupled to 
one environmental oscillator with cji = 1.5wo and A — 0.5a;o at different values of /3. The dashed line indicates the 



entropy in the case when the system is completely thermalised (171 and the dotted line is the maximum entropy 
S'max = lii(2). The Gaussian invariant (50) satisfies the correct properties: initially Axx{to) = 1 and the entropy (44) 
is zero. For zero coupling (A = 0), the two oscillators do not interact and the Gaussian invariant remains conserved, 
leaving zero entropy. For general coupling the Gaussian invariant oscillates between 1 and some value > with an 
angular frequency Cu. The corresponding entropy then oscillates between and the thermal entropy. Only in the limit 
when /3 — > and wq — >■ wi (resonant regime) the maximum entropy ^max = ln(2) is reached (for A^x = 0). 



2. A'' + 1 coupled fermionic oscillators 



In this section we consider the more general case of one system oscillator bilinearly coupled to N environmental 
fermions. In order to find the growth of entropy of the system Eqs. (48) must be solved for the statistical correlators 
with initial conditions ( 49 ) . This can be done numerically. We have used the N = 1 case treated analytically above 
as the test case for our numerical studies. For simplicity the system oscillator couples equally to all the environmental 
oscillators, i.e. Xi = A. If the frequencies of the environmental oscillators are taken in a narrow range away from wq, 
they will effectively behave as a single oscillator, leading to similar plots as figures [l] and [2j 

In figures [3] [6] the system entropy is calculated by taking 50 environmental oscillators with frequencies in the range 
of [0 — 5] X ujQ. The equal couplings to the system oscillator are A^ = A = G.15ajo. In figures |3||4] the environmental 
frequencies are equally spaced, i.e. uii = O.Hujq, i = 1, .., 50, with P — 1 and /3 = 0.1, respectively. The system entropy 
rapidly increases to the value of the entropy in case the system is completely thermalised, see Eq. (17). For higher 
environmental temperature (lower /3) the late time entropy gets closer to the maximum entropy S'max = lii(2), which 
is only reached for /3 — > 0. In general, there will be fluctuations in the entropy due to constructive and destructive 
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Figure 3: System entropy as a function of u)Qt for N = 
50 environmental oscillators. The parameters are: uji = 
O.li XLOo, i = 1, ..,50, A = Q.15tJo and = (a;o)"^ The 
entropy rapidly reaches value of the thermal entropy Sth 
(17 1, indicated by the dashed line. The dotted line is the 
maximum entropy Smax = ln(2). At specific wot there are 
small fluctuations of the system entropy. 
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Figure 4: System entropy as a function of cjot for N = 
50 environmental oscillators. The parameters are: uji = 
O.lixwo, i = 1,..,50, A = 0.15aJo and 13 = 0.1{luo)~\ The 
entropy again rapidly reaches the thermal entropy Sth 1 17 1 
(dashed line), which at such high temperatures (low j3) 
almost coincides with the maximum entropy Smax = ln(2) 
(dotted line). 



interference of the environmental oscillators. Due to the specific distribution of environmental oscillators in Figs. 
[3}|4j the entropy is almost constant with small fluctuations at very regular times when the oscillators interfere in a 
constructive or destructive manner. As the number of oscillators increases these features occur on longer timescales 
and the amplitude of fluctuations decreases. In the limit when — > oo these features disappear altogether. 
In flgures [5||6] the environmental frequencies of the 50 oscillators are randomly selected in the interval [0 — 5] x wo. 
The entropy increases to the same values as in figures |3||4j but due to the random choice of frequencies the late-time 
entropy contains some random fluctuations. As A^ increases, the amplitude of the fluctuations 5n/n decreases. In 
figure [7] we have shown the amplitude of fluctuations in late time statistical particle number for A^ = 25, A^ = 50 
and A^ = 100 and it demonstrates that the size of these fluctuations becomes smaller as A^ becomes larger, and it 
is consistent with the expected behaviour, 5n/n ^ 1/^/N. For infinitely many coupled environmental oscillators the 
statistical particle number, or entropy, at late time becomes more and more stable. 

We have also studied the growth of system's entropy for bilinearly coupled environmental oscillators. For a finite 
number of oscillators the fluctuations in particle number and entropy are large, making it difficult to extract a specific 
growth rate, in particular because in this case the system does not seem to approach exponentially its final statistical 
occupation number. For AA^ 3> 1, when oscillations are small and a relatively stable late time entropy is reached, 
a rough estimate of the growth rate can be made. The rate can be defined as F = n(i)/(rtth ~ ^(i)), ^-nd can be 
estimated as F ^ Awq, plus an additional weak dependence on temperature. The growth rate can be studied more 
accurately in true interacting quantum field theories. In that case one system field mode interacts with infinitely 
many environmental modes through quantum loops, such that the thermalisation of the system field can be quite 
accurately described by the perturbative rate, which in the high temperature regime is propotional to T. At late 
times the system's entropy settles to its thermal equilibrium value. This has been tested numerically for bosons with 
a cubic interaction in Ref. |10j . As with regards to fermions, the quantum field theoretical description of fermionic 
entropy will be discussed in the next section. 

III. ENTROPY GENERATION IN A FERMIONIC QUANTUM FIELD THEORY 

We shall now consider the entropy of a fermionic field theory whose action is given by, 

'S'M = J d'^xL^ , = tp{x)i^'^df_,ip{x) - mtp{x)^p{x) + dnt , -Cint = -j4,{x)^{x) - i){x)j^{x) , (52) 

where ^{x) is a space-time spinor, '4'{x) — ip^^^, 3^jj{x) is a spinorial current, jipix) = j^7° and 7^ are Dirac's matrices 
obeying a Clifford algebra with an anticommutation relation, 

{7^7''} = 2r;^^ ,7^- = diag(l, -1,-1,-1). (53) 
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Figure 5: Same plot as Fig. [5[ but here the environmen- 
tal frequencies have been randomly selected in the same 
interval, i. e. oji £ [0 — 5] x wo • Due to the random choice 
of the environmental frequencies, the entropy randomly 
fluctuates around the thermal value. 



Figure 6: Same plot as Fig. [7| but here the environmental 
frequencies have been randomly selected in the same in- 
terval, i. e. Wi G [0 — 5] X tJo . Again the entropy fluctuates 
around the thermal value much more frequently due to the 
random choice of environmental frequencies. 



(Jn/n 




Figure 7: Amplitude of fluctuations in statistical particle number at late time for N— 25 (purple, dashed), N—50 (red, solid) and 
N—lOO (blue, thick solid) environmental oscillators. Sn/n is defined as 5n/n = (fixx — nth)/nth, where fixx is the statistical 
particle number of the system (45l anduth ~ (exp(/3a;o) + 1)"""^ • For all lines A — 0.1, /3 = 0.5 and the environmental frequencies 
have been chosen randomly in the same range, uji G [0 — 5] x luq . In general, there are fluctuations in the late time entropy due 
to the finite amount of environmental oscillators that constructively and destructively interfere. The plot clearly shows that the 
amplitude of these late time fluctuations decreases as N increases. 



The action (52) is a general Ansatz describing many realistic systems. Examples include: a fermionic field in a heat 
bath of many fermionic degrees of freedom (similar to the quantum mechanical case in section II B I ; one quark flavour 
coupled to other quark flavours through the CKM matrix [31] [55] ; one neutrino flavour coupled to other neutrino 
flavours through the PMNS matrix [23l [24j; but also many systems with true interactions, such as Yukawa type 
jipi ~^ Dij'^ji' with (j) a scalar field. 

In the correlator approach to decoherence [U |5J ISHTO] the Gaussian von Neumann entropy of a system is expressed 
in terms of the Gaussian (statistical) correlators of the degrees of freedom of the system. These correlators are 
derived from the density matrix and they are commonly expressed in terms of those fields that diagonalise the free 
Hamiltonian. The reason is that the time evolution of off-diagonal correlators is zero for non-interacting fields, which 
leads to a simple form of the density matrix. Simple means here that the density matrix can be written as a direct 
product of the density matrices for single degrees of freedom. For interacting theories such as the examples mentioned 
above, the off-diagonal correlators (between the different components of the diagonalised Hamiltonian) are in general 
nonzero and the density matrix has a more complicated form. We will discuss this more thoroughly in the coming 
sections. 

In our trivial example, the quantum mechanical case for free fermions, the Hamiltonian is diagonal because there is 
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only one degree of freedom. Thus the Ansatz ([s]) for the Gaussian density matrix is simple (and already the most 
general one) , and so is the expression for the entropy in terms of the statistical correlator ( 24 ) . For a more complicated 
system such as fermionic fields in 3 + 1 dimensions, which we discuss now, the system is described by a spinor with 
four components for every wavenumber. In order to find the entropy, we first find the fields that diagonalise the free 
part of the Hamiltonian, then use an Ansatz for the density matrix in terms of those fields, and diagonalise it by 
transforming it into an appropriate Fock basis. 



A. Diagonalisation of the Hamiltonian 

As usual, fermions are quantised by employing an anticommutation relation, 

= 6^p5^{x-y), (54) 

(here a, /3 G {1, 2, 3, 4} are spinor indices, which are in other equations suppressed) and thus - just as in the case of 
the fermionic quantum mechanics discussed in section [ll]- they are Grassmannian operators. 

Upon varying the action (52 1 with respect to ii{x) and i^ix), one gets the following operator equations for il^ix) 
and ?/'(a:) , 

ij'^dfj^Tpix) - mip{x) = j^{x) , -id^'iJj{x)Y - mi}j{x) = j^{x) , (55) 

related by hermitian conjugation for real m. The simplest nontrivial case is when the current is generated by a mixing 
mass term. In this case j^{x) — X^fLi iTi'Oj'i'jix) and tpoi^) = "(/"(a:), moo = cf. Eqs. Yipi 



Here we shall study only time dependent problems, and we shall work in a spatial cuBe of volume V, such that it 
is convenient to transform these equations into a spatial momentum space, defined by, 

^(^) = 7f E ^(^' ; V3(fc, ^) = ^ / d'x^ix)e-^^-^ , (56) 

k 

where k = (27r/L)fi, L = V^^^ is the linear size of the cube V, n = (rti, ri2, n.3), rij €E Z and Z is the set of integers. 
With these definitions we then get Eqs. (54) and (|55|) in momentum space. 



and 



{Mk,t),iJ^lik',t)}^6^^pSf:j:,, (57) 
(«7°9t-7-fc-m)V.(fc,t) = j^(fc,t), -idti'ik,t)-/° + ^p{k,t)ik-^-m) ^%{k,t), (58) 



hik,t) - ^ / d'xj^{x)e~-''^-^ %{k,t) = ^ / d^xj^(x)e*'=-" . (59) 



where 



Because the problem at hand is linear, there is no momentum mixing. Since we are interested in time evolution, we 
can work in the helicity eigenbasis, in which 

i,{k,t)=^^h{k,t)®ih{k). j^(fc,t) = ^ih(fc,0®a(fc), (60) 

h=± h=± 

where ^/i(fc) are the two-component helicity eigenspinors, satisfying 

Kh = k- S£,h = Kh , (61) 

where h is the helicity operator in the two-by-two (Bloch) representation of Glifford algebra, which can be defined in 

terms of the helicity operator H as follows, H = k - T. = diag(/i, h) = I2 <E) h, I] — ^^^^^ . By making use of the Bloch 
decomposition of the Clifford algebra in the Weyl/chiral representation, 

7° ^ «) / , 7* ^ ip^ ®a\ 7^ = i-y^^^-f'^^'^ -p^ ® I , (62) 
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where a\ p* {i = 1, 2, 3) are the PauU matrices, equation ([58]) (when multiphed by 7") can be rewritten as 

(idt + hkp"^ - mp^)Mk, t) = p^jhik, t) , (63) 

where tph{k, t) is the two spinor whose components describe the two chirahties and Rh and k = \\k\\. Here we shah 
consider the simpler case when the mass matrix is time independent. ^ In this case a further (orthogonal) rotation, 



R = c^-ip 



R^ c^+ip^s^, R-R^ = I = R^-R, 



TfL Tfl 

tan(20) = — , sin(20) = — , cos(20) 
hk UJ 



hk 



, (64) 



tion ( 63 ) is of the form, 



diagonalises equation (63), where oj = + fn^, C4> = cos(0) and = sin((/)). The resulting (diagonalised) equa- 

{idt + ujp^)i'hik, t) = (c20p^ - S24,p^)Jh{k, t) , 



where 



and we made use of. 



Jh = Rjh = 



jhi 

3h2 



(65) 



(66) 



Rp^R^ 



C20P - S20P 



When Eq. (65) is rewritten in components, we get that the positive and negative frequency modes (particles and 
antiparticlesjobey 



[idt + w)i/'/ii 
{idt - w)V'/i2 



m- hk 

j/il H jh2 

UJ UJ 

hk - rri'- 
— jhi H jh2 ■ 



(67) 



In the absence of currents, the problem is reduced to the diagonal one, and there is no mixing between different states. 
We can define a Fock basis |n^(fc)), which - in the absence of interactions - diagonalises the Hamiltonian. However, 
in general this procedure does not diagonalise the density matrix. Only when the source currents vanish and there 
is no initial entanglement between and "0^2 states the Fock basis simultaneously diagonalises the Hamiltonian 
and density operator. In that case we can define the density operator as a direct product (c/. Eqs. (12) and ( [Tl] )) of 
the density operators for the different fermionic components. In general, however, fermionic interactions (modeled by 



the currents jhi,2) generate mixing between the "0/11,2 and jhi,2 fields, as can be seen from Eq. (67). Therefore, this 
mixing should also be included in the Gaussian density matrix for a fermionic field. 



B. Density operator and fermionic entropy 

Following the previous discussion, a more general Ansatz for the Gaussian density operator for an interacting 
quantum field is 



p{t) = I cxp 1 - ^ ^{k, t)Ah{k, i)*^(fc, t) j , 

V k.h ' 



(68) 



where ^ is defined in (66), Z is the normalisation constant determined by Tr[p] = 1 and 

Ah{k,t) = ( ^^'^^ ) , Xh,ij = Xh,ij{k,t) . 

V ■^h,12 Xh.22 J 



(69) 



When the mass matrix is time dependent, m = m{t) and u) = ui(t), wh ich can occur e.g. during a phase transition in the early Universe, 
then a unitary matrix is needed to diagonaUsc the 2x2 problem (|63|, where now 9 = 9{t) and <j> = <p{t). 
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From now on we suppress the momentum labels and time dependence for the matrix elements Xh,ij- Note that the 
parameters only depend on k = \\k\\ due to the assumed spatial homogeneity and isotropy of the state. In our 
Ansatz (68) the positive and negative helicity states, as well as the different momentum states, do not mix, but the 
"0/11,2 and jhi,2 states do mix through the parameter Ai2- One could consider adding linear terms in ^, in (68) that 
would correspond to a nonzero value of ('I'), (^^). However, if the expectation value is initially zero, it will remain 
zero when the system evolves and Eq. (68 ) is the most general Ansatz for p. The density operator can be diagonalised 
by a unitary transformation U 

Tj _ ( cos 9 -e*'^sin6'\ , . 

\ c-"^sme cosd J ' ^ ' 

with 6 = e{t) and 

cos 20= , ^"'^^^ sin20= , e^* = ^, (71) 

v/(All-A22)2+4|Ai2|2' v/(^ll-^22)2+4|Ai2|2' IA12I' 

such that the diagonalised density operator becomes 

m = lexp \~Y.^''f^{k)Ki{k,i)^'i{k) \ , (72) 



k,h 



where 



#f,(fc) = C/#,(fc,i)= (^ll'+jlJj) , A^,(A:,t) = ;7A,(fc,i)f/^- Y A^)' A,± ^ A,±(fc, i) . (73) 

Here, we dropped the time dependence in "^'^ {k), since these are operators in the Schrodinger picture. The eigenvalues 



are 



A± = \{Xh,ll + A,i,22) ± \\I{Km - Ah.22)2 +4|A/i,i2|^ 



= iTr[Ah] ± \ v/(Tr[A,02-4Det[Ah] . (74) 

In the second line the eigenvalues A/i± are expressed in terms of the Gaussian invariants of the exponent of the density 
matrix, that is, in terms of the trace and determinant of A/j = A/j(fc,t). As in the quantum mechanical case, we 
can identify the statistical particle number Nh±{k) = '0jj-|-(fc)'0h±(fc) and introduce a Fock basis \nh±{k)) (not to be 
confused with the Fock basis |ri^_j_(fc)) that diagonalises the Hamiltonian, discussed above). Of course 

Nh±{k)\nh±{k)) = nh±{k)\nh±{k))) , 7V,,±(fc) = VU±(fc)^'»±(^) • (75) 

The trace of the density operator can now be taken easily, and by demanding Tr[p] = 1 the normalisation is 

Yl (l + e-^"±('=^*)) . (76) 

k,h,± 

Thus the density operator can be written as a direct product, 

p{t)= Y[ph±{k,t), ph±{k,t) = {l-nh±{k,t)) + {2nh±{k,t)~l)Nh±{k), (77) 

k,h,± 

where the average particle number is 

nh±{k,t) = (/^h±(fc)) =Tr[p(t)Nh±(k)] = [1 + exp(Ah± (k, t))]-i . (78) 

The entropy is then (cf. Eq. ([l3|), 

S = ^ Sh±{k,t) , Sh±{k,t) = - (1 - nh±{k,t))\TL(l - nh±{k,t)) - nh±{k,t)\n{nh±{k,t)) . (79) 
kh± 
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The main difference between the quantum fermionic oscillator studied in section |IIB and the (free) quantum field 
theoretic oscillator is that for each fermionic mode there are four distinct states (particles), and hence four contri- 
butions to the entropy for each mode k: two from the two helicity states and two from the particle/antiparticle 
states. Furthermore, Eq. (67) implies that both environmental particles and antiparticles will in general contribute 
to the entropy of the system (anti-)particle. In the non-relativistic limit when fc — > particle-antiparticle mixing is 
absent. This is a property of the particular form of the interaction (52), which for the purpose of this paper we take 



to be an operator valued scalar fermionic current density, which can be generated, for example, by a mass mixing 
term. Other possible interaction Lagrangians that occur in nature include: pseudo-scalar, vector and pseudo-vector 
fermionic currents. For simplicity we shall only consider here the scalar fermionic current. 

Just like the quantum mechanical case the average particle numbers nh±{k,t) (and the entropy) can be expressed in 
terms of statistical correlators. Using Eqs. (68), (74) and (76) one finds 



(i/'LV'/i2) 



d 



1 



1 



h2 



dXhM 
d 

d 

d\hA2 

d 



\nZ ^ -{nh+ + Uh-) + ;:(^^?^+ - nt-) 



InZ ^ -{nh+ + Uh-) - 7;if^h+ - rih-) 



v/(Ah,ii-A,,,22)2+4|A,,a2p 

A/i,ll — \h,22 
V(Vll-A;,,22)2+4|A^,12|2 



InZ = -;j^{nh+ - f^/^-) 



2A 



hA2 



dX 



In^ = ^{nh+ - rih-) 



h.l2 



V22)'+4|A^,12|2 
2A;i_12 

V(Vll-A,,,22)2+4|A^,12P 



(80) 



Here ^phi = iihi{k,t), nh± = nh±{k,t) and Xij = \h^ij{k,t). 
time statistical correlators for the iphi fields 



We can easily relate the correlators above to the equal 



hi 



(81) 



where F^^ij = Fh.ij{k] t; t). The average particle number expressed in terms of the statistical correlators is then 

nh± = ^ - liFh.ii + Fh,22) ± sgn(Fft,22 - Fh m)]^^ [Fh. 22 - Fh,iif + AFh^i2Fh,2i ■ (82) 

The correctness of this expression can be checked when going to the single quantum mechanical fcrmion case, thus 
only keeping for example the V'/ii fields and setting tljh2 fields to zero. For specific h, k there is only one remaining 
particle number, and it agrees with Eq. (22). Moreover, in the absence of interactions and with zero initial mixing, 
the Ansatz for the density matrix (68) becomes diagonal, i.e. A12 = 0. As we stated earlier in the introduction to 
this section, the density matrix then becomes a direct product of different single fermion density matrices. Indeed, 
the non-interacting case gives = i — F^^u, fih- — h — ^/i,22, which agrees with the average particle number for a 



5 — Fh^ii, fih- — \ 

single fermionic oscillator ( 22 ) . Thus the entropy for a non-interacting fermionic field is simply given by the sum of 



the entropies of the components of the diagonalised Hamiltonian, what was to be expected. 

In principle we could have also made an Ansatz for the density operator (68) in terms of rotated fermion fields, for 
example in terms of left- and right-handed fields. Of course, the resulting entropy should not depend on the basis 
in which the Ansatz is made, but one basis may be more convenient than the other. In order to clarify this, we can 
define (just as in the single fermion case) Gaussian invariants of the correlators 



1 - 2n/id 



tanh 



' Xh± \ 



(83) 



The fact that the A^-t are Gaussian invariants of the correlators becomes more clear when introducing a 2 x 2 matrix 
of statistical correlators 



Fh,ii 

Fh,21 



Fh,12 
Fh,22 



(84) 



such that 



A^± = Tr[F^] T v/(Tr[F,,])2 - 4Det[F^] 



(85) 



Both the trace and the determinant are invariant under a change of basis, thus also the expressions for the Gaussian 
invariant are indeed invariant, as are the particle number (82) and the entropy (79). Moreover, because with (74) 
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the eigenvalues X± can be expressed in terms of Gaussian invariants of the density operator, Eq. 
relation between the Gaussian invariants of the correlators and those of of the density matrix. 



(83) presents the 



As a final comment, note that the entropy (79) only gives a limited amount of information. The complete density 



operator (68 1 contains more information and we can separate it into two parts. The "mostly classical" information is 



stored in the spectrum, which we define as 



Spec[/5] = {Xi}: 



(86) 



where Xi (i = 1, ■■,N) are the eigenvalues of —\n{Zp) and can be read off from the diagonalised form of the density 
operator. To be more precise, for an A'^-state system the components of the spectrum are 

A, = (Oi|..(l,|..(OAr|(-ln[Zp])|OAr)..|l,)..|Oi), 

where here p is assumed to be written in diagonal form and the \ni) are the Fock states used to diagonalise the 
density operator. An example of the spectrum for a two-state system can be read off from Eq. (72), where the two 
eigenvalues Xh± (74 1 contain the "mostly classical information". They are related to the averaged particle numbers 
nh± = (e^*** + 1)^^ which is what a late time observer entangled with the Fock states of the system would identify 
with a thermal distribution of fermionic particles. On the other hand, the "mostly quantum" information is stored 
in the off-diagonal components of p, which describe mixing (entanglement) between different states in the original 
(non-diagonal) basis. In section III D| below we show both {Xi} and {fii}. 

The Fock states (75) are a natural candidate for pointer states |TS], which are selected by the environment, and 



in which the system becomes classical, explaining the above term "mostly classical". Hence, these Fock states are 
particularly useful when considering the process of classicalization of a quantum system, and in fact we use them to 
define it. The rate of statistical particle number increase we identify as the rate of classicalization The rotation 
matrix that brings the density operator to a diagonal form has no classical analogue, thereby justifying the name 
"mostly quantum". Of course, even though the Fock states (751 do exist at early times, the system is then not yet 
classical. One therefore needs a more precise definition of when the system becomes classical. 

Zurek states |l8j that pointer states are stable under the influence of the environment, but provides no deeper 
insight into why this is so. We believe that the stability of pointer states can be explained by entropic considerations. 
Namely, in the limit when the number of environmental oscillators A^ becomes large, elements of the reduced density 
matrix in the diagonal Fock basis become stable (up to small statistical fluctuations) because most of the volume of 
the total Hilbert space (of the system -I- environment) corresponds to an almost constant average occupation values 
of the Fock states (75). This represents a quantum generalisation of the ergodic hypothesis. 

As with regards to classicality, there is a notable difference between bosonic and fermionic systems. While, in 
the case of bosonic systems, in the high temperature limit, one can speak of large occupation numbers of certain 
oscillators, yielding a well defined classical field theoretic limit, no classical field theoretic limit exists for fermions 
(simply because fermionic occupation numbers must lie in between and 1). However, that does not mean that there 
is no classical limit for fermionic systems. In the case when the fermionic mass m is very big and the spatial size of 
the fermionic system is large {e.g. when there are many available system states and they are dense in energy), one 
speaks of a classical particle limit, even when the occupation number of each of the states is much less than 1. In this 
case, the number of fermions is well-defined (particle number fluctuations are suppressed) and the Pauli blocking is 
not important. An important example of the classical particle limit is the classical fermionic thermal case, in which 
nic^ ^ fce^, such that fermions get distributed according to the Maxwell distribution, which is of course classical. 
An analogous classical particle limit exists for bosonic systems as well. 



C. Generalisation to N degrees of freedom 

Up to now a density operator p for two degrees of freedom (at fixed helicity h and momentum k) was considered. 



Explicit diagonalisation led to a formula for the "phase space" A's and the entropy (79). The diagonal elements were 
represented by the invariants Tr[F/i] and Det[F?i] of the statistical Greens function matrix. Here we consider a more 
general setting with N degrees of freedom. The (Gaussian) Ansatz for p is now 

^exp{-4>laij^pj) , ^1,..,N) . (87) 



Of course, the classicalization rate is observer dependent and different observers will measure different classicalization rates. For example, 
the position operator (A£^) will perceive a different (typically larger) rate of classicalization. 
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Indeed this is equivalent to a p{9' , 9) in the coherent state representation: 



with 



M,, = (e--^) 



Z 



13 

Detfl + e" 



(89) 



as one can easily see by diagonalising the hermitian (real, symmetric) matrices a, M simultaneously (as in the 2x2 
case considered before). The 2-correlators {'ipl.ipi) are related to the statistical "matrix" F^; and to the aik introduced 
above: 



nu = ii^Ui) - -Ski - Fik{t; t) = -{I-A) 



kl ■ 



(90) 



which is a suitable generalisation to many degrees of freedom of the one degree of freedom result (22). {ip\jpi) can be 
also obtained by differentiating Tr[/5] of Eq. ([87| with respect to —akV- 



d 



da 



kl 



Tr[p] = = Tr 



Tr 



exp(-i/)|ai-,'(/'j) 



1 ' 




J daki 





(91) 



where the first term is just {ip\jpi) and the second term is evaluated as 



-Z 



d 

da, 



kl 



1 d 



Z daki 

a 



Detp- 



da. 



kl 



Trln (l + e"") 



(A; 



Ik 



We have used here the identity Det[A] 

function of the matrix a^-. We thus obtain the generalisation of Eqs. (|22|, 



nki 



exp{tr[ln(A)]}, and that (a/9a„)Tr/(a) 
of Eqs. ^ 

A)ki=[{I + e'^r']^^ 



2^ 



(92) 

(/'(a))^-, where / is some 
(93) 



Diagonalising a, A,F with the same rotation we can use a sum of terms of type Eq. (24 1 for the entropy. This is a 
trace and rotating back inside the trace we obtain 



S 



-Tr 



A. 



In 



In 



(^) 



-Tr [(I — n) ln(I — n) -|- n In n] 



(94) 



For this formula to work the eigenvalues of h must lie in the interval [0,1], which is indeed the case for fermionic 
systems. The result (94 1 we have also obtained using the replica trick for the density operator p in the coherent 
representation (88 1 without a diagonalisation procedure. Appendix [P] contains the calculational details. Of course in 
order to evaluate (94), diagonalisation of A, a is again the fastest method to obtain the entropy. In the next section we 
discuss the growth of entropy for Dirac fermions mixing through a mass matrix, where we will also show the spectrum 
oi p (86) as well as the average particle number. 



D. Fermion mass mixing 



A simple model for interacting fermions is a model of different fermion species mixing through a mass matrix. 



Similar to the quantum mechanical action of bilinearly coupled fermions (41 ), the action for fermion mass mixing 



(95) 
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with 



N 



1=1 



i=l 



(96) 



From now on we assume that the mass mixing parameters are real, mg^ — rriiQ. Next, we foUow the same steps as in 
section III A we first transform all fields to momentum space as in ( 56 ) , then go to the helicity eigenbasis using ( 60 1 , 



and finally rotate the fields as in Eq. (66) (with a different rotation matrix Ri for different species). The resulting 
action is 

^[4,^5 J] - f dtJ2\'^l.hik,t) {idt+ujap'') 'i',Mik,t)+J2'^l.h^^^^^ {idt+u;^p^) ^^../.(fc, t) 



N 



(97) 



where wq = y + II^IP ^^"^ ~ \l +11^11^- ^ — ^(^) is the rotation matrix that diagonalises the free part of 
the action for ip^, whereas Ri = R{di) diagonalises the action for -0^. . When the masses of different species are the 
same, mo = rrii, the rotation matrices will be the same, R = Ri. In general the interaction term is 



Rp^Rj = R^p^R^ = cos(6i + ei)p^ - sin(6i + ei)p^ , 



with the 0, 9i defined as in Eq. (64 1. The equations of motions follow directly from the action (97) 



(98) 



N 



{idt + ujQP^)^.c,h{k, t) = Y^ mQ,Rp^Rj^q^^h{k, t) 



i=l 



[idt + w,p^)*g. ji(fc, t) ^ mQiR,p^R^^^^h{k, t) 



N 



+^o^L,(fc,i)p') = Y.'^^^^Lhik.t)R,p^R' 



1 uT 



i-^dt^l^ik. t) + uj.^l uik. t)p') = mo,VU(fc, i)Rp 



(99) 



Note that each line consists of two equations for the two components of the spinors '^x,h and Vl/q. h- Differential 
equations for the statistical equal-time correlators can be derived from the equations of motions (99). It is only 



necessary to derive the statistical correlators for fields with the same helicity. Remember that there is no helicity 
mixing in the action. For the system alone, the statistical correlators obey, 



N 



idtFxx,hn = X! [cos(6' + 9i) {Fq^^Mi - F^q^^hii) - sin(6' + 61.,) {Fq.x^hii - F^qijin)] 

i=l 
JV 

ldtFxx,h22 = ^ moi [C0s{e + 9i) {Fq^x.hl2 - Fxq^,h2l) + Sm{0 + e.i) {Fq.x,h22 " ^'xqi,ft22)] 
i=l 
JV 

{tdt + 2ujo)Fxxjii2 = ^ mai [cos{9 + 9i) {Fq^x,h22 ~ Fxq^,hii) - sin(6' + 6*,) {Fq.xM^ + F^q^jiu)] 



(100) 



Here we have used a shorthand notation, with 



1 



Fxxjiran = Fxx,hm7i{k]t]t) = -^^iii^ x Jirn{k , t) , ^,^^{k , t)]) , 



Tn,n = 1,2, 



etc . 



(101) 
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Note that F*^ 



xx,h.nm, ^'xqi.hTnn = Fqixjinm and F*,^, ,^^^ = Fq-q.^hnrn- With these relations the remaining 
equation for Fxxji2i in (100) can be found easily by taking the complex conjugate. Similarly, the equations of motion 
for the environmental correlators are: 



[idt + {Uji - Ulj))Fg^q-^hll 
[idt - {cJi - UJj))Fq.q^^h22 
{idt + [uJi + UJj))Fq^q^^hl2 



rriQ, [cofi{e + ei)F.xq^^h2i - s\n{e + 6,)Fxq^^hii\ 

- niQj [008(6* + 0-j)Fq^x^hi2 - sin(6' + Oj)Fq^xjiii] 
moi [cos{0 + 9i)Fxq^^hi2 + sin{9 + e,)Fxq.j,22] 

- moj [cos{0 + Oj)Fq^x^h2i + sm{e + 9j)Fq.xj,22] 

TUq, [C0s{e + 9i)Fxq^M22 - Sin{9 + ei)Fxq.Jil2] 

- moj [cos{9 + 9j)Fq^,x^hii + sin(6' + 9j)Fq.x^hi2] 



(102) 



Again, the remaining equation for Fq.q., h2i can be obtained by complex conjugation of the third line above. Finally, 
the system-environment correlators obey the equations 



{idt + (wo - ^i))Fxq^Mll = 



{idt - (wo - UJi))Fxq^ 



h22 



{idt + (wO + ^^i))Fxq„h 
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{idt - (wo + UJt))Fxq^^h21 



N 

^ moj [C0s{e + 9j)Fq^q.^h21 - Sm{9 + Oj)Fq.q^^h.ll] 
3 = 1 

- moi [cos{e + 9,)FxxM2 - sm{9 + 9,)Fxx^h.ii] 

N 

^ mOj [C0s{e + 9j)Fq^q. Jil2 + Sin{9 + 9j)Fq.q^,h22] 

3 = 1 

- [cos{e + 9^)Fxx,h2i + sin(6l + 9{)Fxxj,22] 

N 

^ moj [C0S(6' + 9j)Fq^q.^h22 - Sin(6' + ej)Fq.q^^hl2] 

i=i 

- moi [cos{e + 9i)Fxx^hii + sm{9 + 9i)Fxxjii2] 

N 

^ mOj [C0s{e + 9j)Fq-q.^hll + Sin{9 + ej)Fq.q^,h2l] 
3 = 1 

- moi [cos{9 + 9i)Fa:xM22 - sm{9 + 9i)Fxx^h2i] ■ 



(103) 



Taking the complex conjugate of these equations gives the final equations of motion for the environment-system 
correlators. This results in a closed system of {N + 1)^ x 2^ x 2 equations for the correlators of the components 
of {N + 1) coupled 2-spinors at different helicities. These coupled first order differential equations can be solved 
(numerically) with initial conditions corresponding to environmental oscillators in chemical equilibrium: 



Fxx,hiii^'^ ^o; to) — FxxM22{k; to; to) — 



1 



,hii{k;to;to) = % - tanh 



'^q,qj,h22{k;to;to) = (5ij-tanh 



l3{uj^ - ^i) 



l3{uj^ + ^i) 



(104) 



and all others are initially equal to zero. According to Eqs. ( 104 \ at t — to there is no mixing between the different 
components of the 2-spinors in the helicity eigenbasis. Remember that these '1,2' components are the fields that 
diagonalise the Hamiltonian in a non-interacting theory; they are the positive and negative frequency states, or 
particles and antiparticles. Note that the initial state ( 104 ) allows for nonvanishing chemical potentials fii for the 
environmental fields. The chemical potentials have an opposite sign for particles and antiparticles. Moreover, we have 
assumed that initially there is no mixing between the system and the environment. The physical picture is therefore 
that initially there is no mass-mixing between the different fermion species, but at t = to the coupling is switched on 
and the entropy of the system can grow. That is, we consider only the entropy of the system 



Sx{t) - J2 

{k,t) 

kh± 

= E 



kh± 



1 + ^xx,h±{k, t) ^^(1 + A. 

xx.h:t 

{k,t) 



1- A 



xx.hzt [ 



{k,t) ^ f 1- Axxj,±{k,t) 



In 



(105) 
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Figure 8: System entropy as a function of mot for N — 1 
environmental field and zero chemical potential, /^i = 0. 
The parameters are mi = l.lmo, moi = 0.5mo and /3 = 
(mo)~^. The entropy is expressed per (^)^- 




Figure 9: System entropy as a function of mot for N = 1 
environmental field and zero chemical potential, fii — 0. 
The parameters are the same as those m Ftg. [£[ hut the 
temperature is twice higher, j3 — 0.5(mo)^^. The late 
time entropy is approximately 8 times higher than for j3 = 
(mo)~^, which supports the scaling of late time entropy as 



where S2;2:,/i±(||fc||, i) is the system entropy per fermionic degree of freedom, i.e. for a state with quantum numbers 
k,h,±. The Gaussian invariants Axx,h± £^re those defined in Eq. (85), with the subscript xx indicating that only the 



system correlators are used. Due to the loss of information (assuming environmental correlations are inaccessible) the 
system decoheres, leading to an increase in entropy for the system. 

In the case of = 1 environmental fields we have numerically solved the 16x2 equations for the statistical correlators 
of particles and antiparticles of the system and environment at different helicities. In the numerical procedure a smooth 



selection of modes k has been made, separated into spherical bins of size Afc. The total system entropy ( 105 1 is then 
calculated as 



Sxit) = ^ Sxx,h±{k,t) 



kh± 



fc/Afc=0'i=± ± 



k 

mo 



Afc 

mo 



h±{k,t) 



where V is the volume of the system and Sxx,h± {k, t) is the (average) entropy per degree of freedom in a spherical 
bin with a momentum ||fc|| and a width Afc ^ k. The maximum mode has been chosen such that /3fcmax S> 1, since 
the inclusion of higher modes does not significantly change the total entropy. In Figs. |8] and [9] the entropy density 
for the system (in units of the inverse Compton wavelength cubed, — (mo/(27r))^ has been plotted for zero 
chemical potential, same mass and mass-mixing parameters, but different temperatures. In the absence of a chemical 



potential the particles and antiparticles evolve completely separately, i.e. Fx 



,hl2 



and Fxx h2i are zero. Because the 



initial conditions are identical the statistical correlators Fxx,hii and Fxx,h22 behave equally, and so do the Gaussian 
invariants Axx,h± of Eq. ( |85[). Thus the total entropy is simply four times the entropy calculated from a single 
Gaussian invariant. In Figs. |8| and |9] the total system entropy increases due to interactions with the environment. 
After some time it fluctuates around an equilibrium value. This late time entropy should scale as j3^^ = {k-QTY in 
the relativistic limit where k^T/niQ ~ l/(/3mo) ^ 1. Comparing Figs, [s] and [o] this appears to be the case. 
The fluctuations in late time entropy in Figs. |8] and [9] are rather large. The reason is that the bilinear coupling is 
not a true interaction term: each system fleld mode is only coupled to A'^ = 1 environmental mode. Due to the 
unitary evolution, energy flows back and forth from the system to the environmental oscillator, resulting in large 
amplitude oscillations. As mentioned at the end of Sec. II B 2[ in the case of a true interaction each system mode 
couples effectively to infinitely many environmental modes, leading to an efficient thermalization of the system with 
an expected rate that is to a good approximation given by the perturbative rate, just as in the case of bosonic field 
theory [TO]. 

Next, a distinction can be made between particles and ant iparticles by introducing a nonzero chemical potential 
/ii for the environmental fermion species ipq^ 



In Figs. 



10 



13 



the particle/antiparticle number densities (82) and the 
spectrum of the density operator (|86[) are shown at different values of /3. We have taken here the case of a system 



field interacting with 10 environmental fields, which have masses distributed around the system mass mo. In general, 
the (anti)particle number oscillates between the initial value and the value for perfect thermalisation, when the 
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nh± 
0.5 r- 




Figure 10: Average (anti)particle number nh± for the 
mode k = mo as a function of mot for N = 10 environ- 
mental fields. The environmental masses are distributed 
as mi = (0.5i — 0.25) x mo, or mi = 0.25mo,m2 — 
0.75mo, ..,mio = 4.75mo, and the couplings are all equal 
moi = 0.2mo. The inverse temperature is 13 = (mo)""*^ 
and the (equal) chemical potentials are Hi = mo. Due 
to the chemical potential the antiparticle number density 
fih- (solid blue) is suppressed with respect to the parti- 
cle number density nh+ (solid red, thick). The dashed 
lines indicate the (anti)particle number densities for per- 
fect thermalisation, (nth)h± = (e^'"^''^'+l)~"'^. Note that 
there is no distinction in particle number for + and — he- 
licity states because helicity mixing is absent. The dotted 
black line is the maximum fermionic particle number in 
the limit when /? — )■ 0, nh± 0.5. 



Figure IT. Average (anti)particle number nh± for the 
mode k = mo as a function of mot for TV = 10 envi- 
ronmental^ fields. The parameters are the same as those 
but the temperature is higher, (3 = 0.5(mo)^"'^. 
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m Fig. 

Both the particle (solid red, thick) and antipartic le ( solid 



blue) number densities are larger than in Fig. 10 but 
the relative increase of the antiparticle number density is 
bigger. 




12: Spectrum of the density operator for the mode 
k = mo as a function of mot for N = IQ environmental 
fields. The spectrum, defined in Eq. (861, are the \h± 
of E q. ( 83 1 . The parameters are the same as those m 
Fig. 10 Due to the chemical potential \h- (solid blue) is 



greater than \h+ (solid red, thick). When the system is 
completely thermalised the spectrum is Xh± = /3(aJoTMi); 
indicated by the red and blue dashed lines for particles 
and antiparticles, respectively. Note that the initial value 
of Afi± IS infinite as the initial particle number is zero. 



Figure 13: Spectrum of the density operator for the mode 
k = mo as a function of mot for N — 10 environmental 
fields. The parameters are the same as those in Fig. | jg[ 
but the temperature is higher, (3 = 0.5(mo)~"'^. As tem- 
perature increases, the difference between Xh+ and A;,_ 
becomes smaller. 



system fermions have the same temperature as the environmental fermions, approximately the initial temperature of 
the environment. Moreover, Figs. [TO] and [TT] clearly show that, for positive particle environmental chemical potentials. 
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the system antiparticle number density is suppressed with respect to the particle number density due to the nonzero 
chemical potentials. As usual, for higher temperatures (lower /3) the particle numbers are closer to the maximum 
fermionic particle number fimax for /3 -> 0, but the relative increase of the antiparticle number with respect to particle 
number is greater. 



The spectrum of the density operator (86 1 is shown in Figs. 12 and 13 and is related to the particle number as in 
Eq. (75). The eigenvalues Xh± of the exponent of the density operator are initially infinite (corresponding to zero 
(anti) particle number), but oscillate on top of its thermal value of Xh± = l3{ujo T Mi) later time. When more 
environmental fields are added, the oscillations are damped and the eigenvalues, and thus the (anti)particle numbers 
move closer to a constant. This is similar to what happened in the quantum mechanical case, see Figs. |3}|6j 

IV. DISCUSSION 

In this work we provide a quantitative description of the entropy of quantum mechanical and quantum field theoretic 
fermionic systems, which here consists of one system oscillator (or field) and N environmental oscillators (or fields). 
In our correlator approach to decoherence the observer is assumed to be sensitive only to the 2-point correlators of 
the fermionic system oscillator (or field). In that case the reduced density operator of the system is Gaussian and 
the corresponding Gaussian entropy can be explicitly calculated in terms of the correlators. We have done this for a 



one-dimensional fermionic harmonic oscillator (24), for a fermionic quantum field theory (79) and for the general case 
of N fermionic degrees of freedom (94). We have demonstrated that the Gaussian density matrix singles out a Fock 
basis in which it becomes diagonal. The Fock basis defines statistical particle number, whose dynamics can be used 
to define how a system evolves from quantum to classical. In that sense this Fock basis defines pointer states. An 
observer which measures statistical particle number, when it gets entangled with these pointer states, will experience 
complete decoherence. The decoherence rate associated with that observer can be defined to be the classicalization 



rate. In Sec. II B we make a rough estimate of this rate for the quantum mechanical case studied here, but we leave 
a more detailed study of the classicalization rate in realistic fermionic quantum field theories for future work. 

For simplicity in this work we have considered fermionic problems where the system couples bilinearly to the envi- 
ronment. This problem has the advantage that it can be solved exactly by numerical methods. We have demonstrated 
that the dynamics in general leads to an increase in entropy of the system. When the system couples strongly to 
the environment in a thermal state at temperature T, at late times the system's entropy reaches its thermal value at 
the same temperature T. Furthermore, in the field theoretic case we have shown that, when environmental fermionic 
fields are in a chemical equilibrium with the common particle chemical potential, the system field will eventually reach 
chemical equilibrium with the environment. 

While in this work we focus our attention on the study of exactly soluble Gaussian systems with bilinear couplings, 
fermionic systems occuring in Nature are usually not of that type. An important example is a relativistic quantum 
field theory with Yukawa interactions ([2|, where the scalar field is the Standard Model Higgs, a candidate for which 
has recently been discovered [28l [29], or the inflaton in inflationary models. A more sophisticated treatment of the 
scalar field in the Yukawa interactions is desirable, and one can foresee solving nonlinear, perturbative, Kadanoff- 
Baym equations for the fermionic and scalar fields, whereby scalar thermal fluctuations are also taken into account. 
Analogous equations have already been tackled within a certain approximation scheme for the bosonic case in [Tll^llO). 
We intend to address the analogous problem for fermions in future work. 

Another interesting extension of this work would be to study the effects of CP violation by adding coupling to a 
pseudo-scalar fermionic current with time (or space) dependent mass mixing terms (thus emulating phase transitions 
in the early Universe). In this case, the CP violation would induce a difference between the particle and antiparticle 
numbers, which in the massless limit becomes the axial vector current. Studying how this axial current depends 
on the environmental temperature in the presence of a non-adiabatically changing mass would allow for a better 
understanding of baryogenesis and leptogenesis sources 
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Appendix A: Bosonic density operator and entropy 

In this appendix we calculate the invariant (phase space) area and entropy for a quantum mechanical system of 
one bosonic degree of freedom with position operator (j) and momentum operator tt, based on the approach in Ref. 
0] . The Ansatz for the bosonic density operator is 



PB[t) = -exp 



1 



(Al) 



where {., .} is the anticommutator and ck,/3,7 are real time dependent parameters. By defining bosonic creation and 
annihilation operators 



Ob 



[0- 








) + l — TT 






a 





1 " 




(j) — l — TT 






(7 





the density operator can be written in diagonalised form 

PB(t) = ^exp(-(TiVB), 



Z' = 



z 



(A2) 



(A3) 



The bosonic particle number is defined in the usual way N-^ — ajjCiB, with a corresponding Fock basis Inn) defined 
through A^bI't-b) — n-Q\n-Q). Using this basis to take the trace and demanding that Tr[/5B] = 1 we find 



Z' = Tr[exp(~aiVB)] = ^ (tibI exp(-fT7VB)|nB) = ^ 



nB=0 



TlB=0 



1-e- 



(A4) 



The average particle number is 



(TVb) = Tr[pBA^B] 



riB 



(A5) 



which indeed agrees with the Bose-Einstein distribution for a thermal state if we identify a = E/{kBT), where is 
the Stefan-Boltzmann constant. The Gaussian correlators are obtained from /5b as 



d 

71-2) ^ _2— InZ = (fiB 



da 



1X7 



27 cr 



1\ a 



= nB 



-({0,7r}) = (nB + 



2J a 



2/ cr 



The statistical correlator for bosons is 



(A6) 



F4t;t') = -mt)A{t')}}, 



(A7) 



which we use to define a Gaussian invariant 

,1 



A^(t) =4 



2)(7r2)-(-{0,vr})2 



= 4 [F^{t;t')dtdt'F^{t;t') ~ {dtF^{t;t')f 



(A8) 



A0/2 is the phase space area occupied by a Gaussian state in units of h [7 . In a free theory — 1 and conserved, 
whereas it increases for interacting theories. Using the correlators (A6) we find 



A^(i) = l + 2nB(t) 



tanh(cr/2) ' 



(A9) 



which presents a relation between the invariant phase space area and the Gaussian invariant of the density matrix cr, 
and should be compared to the expression for fermions ( 23 ) . Finally, the bosonic entropy is 



S4, = -Tr[ln/5B] 



1 + A^ ^^ l + A^ 1-A^ ^^^ l-A^ 
2 2 2 2 



In — ^^i^ — — - In — - = (1 + ^b) ln(l + tt-b) — "-B InriB 



(AlO) 
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Appendix B: Fermionic shift and diagonalisation 



The lagrangian for an interacting fcrmionic oscillator is given by (see Eqs. (24) and (39)), 

= tp\tdt - uj{t))i} - jlip - ^'^j^, . 

This implies the equations of motion, 

{idt - = , {-idt - u}{t))ip^ = . 

One can easily construct the free field solution (in the absence of currents) , 



in terms of which we can express the (free) retarded and advanced Green functions ( 20 ) as, 

iS'„{t;t')^^ie{t-t')c-'I>'^^ , iS^[t-t') = id{t' -t)e'i>'^^ 



With a help of iSl^ we can solve the general fermionic operator equations (B2), 



Jo Jo 
Now, for free fields '!/'o(0 — e^^-'o "'''^ 7/10(0) and the related ipoit) "^"^ have, 

Noit) = i^lWMt) = i^l{o)Mo) = ^o(o) 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



which implies that in this (source-free) case the density matrix (11) does not evolve in time, such that a = const., 

fi = const, and also the entropy (13) S — const. 

In the case when there is a nonvanishing current source, j^(t) 7^ 0, one can think of the full solutions (B5) as a 
suitably shifted "00 (i) (analogous to the bosonic Glauber's coherent states), and therefore one can write the density 
operator in terms of the shifted fields 



(B7) 



Eq. ( B6 ) then implies that for that density operator the (von Neumann) entropy is conserved, as it should be. 



Appendix C: Exact entropy for two coupled fermions 



For one environmental oscillator the equations of motion for the statistical correlators Eqs. (48) become 

idtF^^{t;t) = --XAF{t;t) 

idtAF{t;t) = {uo - L^i)F+{t;t) - 2X{F,,{t;t) - Fg,it;t)) 
idtF+{t-t) = (wo-wi)A^^(i;i) 
idtFqq{t-t) = XAF{t-t) , 



(CI) 



where AF(t; t) — Fxq{t; t) — Fqx{t; t) and F^{t; t) — Fxq{t; t) + Fqx{t; t). The statistical correlator can be solved from 
the first line 



Fxx 

(t; t) — Fxx{to: to) 



A 



f dt'AF{t';t') = --^ f dt'AF{t': 

J to ^ Jto 



t'), 



where we have used the initial conditions ( |49| . After acting with idt on the second line of (CI) one finds 

(92 + (wo - ^1)' + 4A2)Ai^(t; 0=0. 



(C2) 



(C3) 
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This equation can be solved with initial conditions for the correlator AF{t;t) itself from Eqs. (491 and for the first 
time derivative of the correlator AF{t;t) from the second line of (CI). This gives 



2A 

AF(t; t) = I— HE sin[u){t - to) 

CO 



(C4) 



where 



UJ = ^/{uo - wi)2 +4A2 
1 



TIE = 



1 1 

- = - — - tanh 
12 2 



Inserting the solution ( C4 ) in ( C2 ) one obtains the Gaussian invariant from Eq. ( 45 1 , 

2 



A^x{t) = l-2nE 



2AV 



sm 

UJ J 



(C5) 



(C6) 



Thus, the entropy of the system can be analytically calculated using (44). 



Appendix D: Entropy via the replica trick in coherent state basis 



Here we calculate the Gaussian von Neumann entropy ([T]) using the density operator in the coherent state basis 
(33). For convenience we use the exponentiated form of the elements of the density operator (see Eq. (29 1), 

p{e',e;t) = {9'\p\0) ^^exp{e'M0) , (Dl) 

where M = = e^" and Z = = 1 + e^" = 1 + M. By making use of the replica trick the entropy can be 
expressed as 



S — — Tr[/51n(/3)] = — lim 



Tr[p"+i - p] 



(D2) 



The trace is defined in (30). By inserting n unity operators (31) in (D2) and using {9'\9) = exp{6'6) and Eq. (Dl) 
one finds 



Tr[p"+i] = J dOdOexpiee) f[ J 



Z 



= Z 



-ri-l 



= 1 

dOdd exp 
-"-i(l + A/"+i) 



exp( 



n 



de^'^de^'^ exp(-^(*)6i('); 



exp(6iM6i(i)) X exp(6l(i)A/6l(2)) x .. x exp(6'(")Af6') 

(D3) 



The Grassmann integrations have been performed explicitly in the last step. The resulting entropy (D2) becomes 

1 + A'/"+i 



S — — lim 



o\{l + Af 



M 
1 + M 



In Af + ln(l + Af) = -(1 - n) ln(l - n) - nlni 



(D4) 



which is indeed the same as the entropy derived earlier using the Fock basis ([T3|). 



The previous derivation can be generalised for N fermionic degrees of freedom. In that case (see Eq. (88)) 



with 



(D5) 



Z = Det[I + e""] = Det[I + Af] 



(D6) 
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Again the replica trick (D2| is used to obtain the entropy. The trace is of course taken over all N fermionic degrees 
of freedom, which are also included in the unit operation. After inserting the unity operators in ( D2 ) one finds 

N 



Tr[p"+i] = n / d9ad9a exp (^^e,) x f[ f S'^Uo^^^ exp (-e\^^e^^- 

a exp [OaOa 



a=l 

X 



N r „ 

Y\ / d^(")(i6'(") exp (-^<;"^6l(")) 



n[/ 



n /de^d^Wexp(-0-W0«) x...xn /rfe'rfe)exp(-e^e 
b=i -I c=i L"' 



x exp \Y,hM,/p X exp 



fc,/ 



(D7) 



To avoid confusion the summations have been written out explicitly. In order to perform the Grassmann integrations 
the following identities prove to be useful, 



Det[l + M] = ]^ / d0ad0a expiMa) exp{J2 O^M^Je,) 

N 



id 6=1 '-•^ 



rfe-Wde«exp(-0-«eW) 



exp(^^,M,fe0W) X expC^el^hiijej) 

I, J 



i,k 



Applying these to Eq. (D7| we find 



Tr[/5"+i] = Dot [I + A/"+i] 



The entropy (D2|) becomes 
jS* = — lim 



Dct[I + M"+i] 



"^0 1 (Dct[I + Af]) 



n+l 



- 1 S = Tr 



M 

'l + M 



lnM + ln(I + M) 



(D8) 
(D9) 

(DIO) 

= Tr[-(l-n)ln(l-fi)-nlnn] , (Dll) 



which is indeed the entropy derived in Eq. (94 1. Thus, by using the density operator in the coherent state basis in 
combination with the replica trick, no diagonalisation of the density operator is required in order to find the entropy. 
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